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Abstract. Interlayer tunneling in graphite mesa-type structures is stud- 
ied at a strong in-plane magnetic field H up to 55 T and low tempera- 
ture T = 1.4 K. The tunneling spectrum dl/dV vs. V has a pronounced 
peak at a finite voltage Vo. The peak position Vb increases linearly 
with H. To explain the experiment, we develop a theoretical model of 
graphite in the crossed electric E and magnetic H fields. When the 
fields satisfy the resonant condition E = vH, where v is the velocity 
of the two-dimensional Dirac electrons in graphene, the wave functions 
delocalize and give rise to the peak in the tunneling spectrum observed 
in the experiment. 
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1 Introduction 



Over the last decade, interlayer tunneling spectroscopy was developed to measure 
energy gaps in high-temperature superconductors and charge-density-wave materi- 
als [112] , It was also used to study magnetic-field- induced charge-density waves in 
NbSe3 [3] and graphite A . This latter effect has orbital origin and generally exists 
for a magnetic field oriented perpendicularly to the conducting layers. In contrast, 
the interlayer tunneling in a parallel magnetic field can be utilized to obtain infor- 
mation about the in-plane energy spectrum of the carriers in the layers. Tunneling 
in a GaAs/AlGaAs bilayer subjected to parallel magnetic and perpendicular electric 
fields was studied in [516] . Tunneling experiments between two graphene sheets in the 
graphene/insulator/graphene heterostructure were recently reported in [7J. The effect 
of the in-plane magnetic field was studied theoretically for graphene multilayers [5] 
and for a thin film of a topological insulator theoretically |9I10| and experimentally 
[llj . In this paper, we present an experimental and theoretical study of graphite in 
strong in-plane magnetic and out-of-plane electric fields. 



Will be inserted by the editor 

-V _^±v 




Fig. 1. Schematic view of the mesa-type structure. 
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Fig. 2. Interlayer tunneling spectra of the graphite mesa-structure at various in-plane mag- 
netic fields at T = 1.4K. 



100 

80 
60 
40 
20 

i ° 

—- -20 
> -40 
-60 
-80 
-100 




5 10 15 20 25 30 35 40 45 50 55 60 

H CO 



Fig. 3. Dependence of the voltage Vo of the tunneling conductance peak on the in-plane 
magnetic field H . 



2 Results 

Mesa-type structures were fabricated by double-sided etching of thin graphite flakes 
using focused ion beam [TJ]. We used high-quality natural graphite of NGS Natur- 
graphit GmbH Co. The mesa is shown schematically in Fig. [T] Electric current flows 
from one part of the crystal into another across the layers in the region of the mesa. 
The mesa is typically of 1 micron size and contains a few tens of graphene layers. 
Because of a very high interlayer conductivity anisotropy (<r±/cr\\ ~ 10 4 at low tem- 
peratures), the applied voltage V drops mostly on the mesa. The experiment was 
performed in pulsed magnetic fields at the National Laboratory for High Magnetic 
Fields in Toulouse. We developed a system of fast data acquisition, which allowed 
us to collect about 10 3 I-V characteristics, each containing 10 3 points within each 
magnetic field pulse of 400 ms duration (3J. Figure [5] shows a set of dl/dV spectra vs. 
V for various in-planc magnetic fields. For magnetic fields greater than 20 T, local 
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Fig. 4. Schematic view of two graphene layers and the directions of the applied electric and 
magnetic fields. 

maxima start to develop for both polarities of the bias voltage symmetrically. Both 
the height and the voltage Vo of the peaks grow with the magnetic field H. Figure [3] 
shows that the dependence of Vq vs. H is close to linear. 



3 Discussion 

In this section, we consider a model of graphite in crossed electric and magnetic fields 
in the geometry shown in Fig.0] The magnetic field H = Hy is parallel to the layers, 
while the electric field E = Ez is perpendicular, and the distance between the layers 
is d. We shall demonstrate that, when the resonance condition for the fields 

E = vH (1) 

is achieved, the electronic wave functions become delocalized in the z direction and 
give rise to a peak in the interlayer differential conductance, as observed experimen- 
tally in Figs. [2] and [21 The parameter v = 10 6 m/s in Eq. ([XJ) is the velocity of the 
two-dimensional Dirac electrons in graphene. 

The Lorentz force from the in-plane magnetic field H acting on the interlayer 
electron tunneling shifts the in-plane momentum p x of the electrons on the n-th 
graphene layer to p x — qn, where q — eHd [8|10j . The uniform out-of-plane electric 
field E changes the potential energy on the n-th layer by un, where u = eEd. Thus, 
using the approach of Ref. [8], the Schrodinger equation for graphite in the crossed 
electric and magnetic fields can be written as 

va x (p x - qn)W n + tI A {<P n _ 1 + W n+l ) = (e - un)W n . (2) 

Here, the Pauli matrix a x acts on the spinor wave function <?Vi = which 
has components on the A and B sublattices of the n-th graphene layer. The matrix 
I A = (l+er z )/2 and the amplitude t = 0.4 eV describe the interlayer coupling between 
the carbon atoms, which lie on top of each other in the Bernal-stacked graphite lattice. 
For simplicity, we set the in-plane momentum component p y = below; however all 
calculations can be easily generalized for p y ^ 0. In the limit of zero interlayer coupling 
t — > 0, the energy spectrum is given by a series of Dirac cones 

e m = ±v(p x - qm) + um, (3) 

which correspond to the wave functions localized on m-th graphene layer. Evolution 
of the spectrum ([3]) with the increase of the electric field is qualitatively illustrated 
in Fig. [5] Each Dirac cone consists of two branches corresponding to the left-moving 
and right-moving electrons, referred to as the L-movers and R-movers in the rest of 
the paper and labeled in Fig. 02a) . For u = 0, the Dirac cones are shifted horizontally 
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(a) u = 



P, 




(b) u<vq 




Px 



(c) u = \q 



Fig. 5. Schematic energy spectra of graphite in crossed magnetic and electric fields (rep- 
resented by q and u) in the limit t — > at different values of u. (a) The Dirac cones have 
the same energy for zero electric field at u — 0. (b) With an increase of u, the Dirac cones 
shift vertically, (c) When the resonant condition u = vq is met, the R-branches of the Dirac 
cones align and become degenerate. 



by q, as shown in Fig. [S[a). For u ^ 0, the Dirac cones are also shifted vertically, as 
shown in Fig. [SJb). When the fields E and H satisfy the resonant condition u = qv, 
which is equivalent to Eq. (fl}, the right- moving branches of the spectra align and 
become degenerate, as shown in Fig.|5[c). 

For a more detailed analysis, let us perform a unitary transformation W n = 
e -io- S 7T/4 ! p7 tQ tlle bagis Q £ t j ie L an( j moverSj \p' n — ipLy Then, Eq. ^ becomes 

va z (p x - qn)V' n + 1(1 - a*)^'^ + = (e - un) S£. (4) 

If we temporary neglect the matrix a x , then Eq. (HJ decouples for the R and L modes: 

[vp x + n(u - ^ + -(^_! + </£fi) = (5) 
[»p x + n(u + qv)} ift + 1 + - e^. (6) 

Solutions of Eqs. ((SJ and (O can be enumerated by an integer to [5], so that the 
eigenvalues are given by Eq. ([3]), and the eigenfunctions are 

^ L = Jn-m ( — —) , (7) 

where J m (x) is the Bessel function. The spectrum © is given by the two sets of curves 
= vp x + m(qv — u) and = — vp x + m(qv + u), which have linear dispersion in 
p x . Thus, a combination of the electric and Lorentz forces shifts the linear dispersion 
by (u — qv) for the R-movers and by (u + qv) for the L-movers. This results not only 
in different spacings between the R and L modes in Eq. but also in different 
localization properties of the corresponding wave functions. The wave functions (JZJ) 
are localized in the z direction on a finite number of layers of the order of tj\u^qv\ . So, 
with an increase of the electric field, the L (R) movers become more (less) localized. At 
the resonance condition u = qv, the R-movers experience zero net force, so their wave 
functions become delocalized, and the spectrum acquires a tight-binding dispersion 

e R = vpx + tcos(k z ), (8) 

where the out-of-plane momentum k z is a good quantum number. 
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Fig. 6. Numerically calculated energy spectra and wave functions for a finite system of 150 
layers. The top panels show the energy spectra, whereas the bottom panels show the wave 
functions corresponding to the states marked by the (*) and (x) symbols in the top panels. 
The left and right columns correspond to the non-resonant 2u = qv and resonant u — qv 
cases, respectively. 



The numerically calculated energy spectra of the full Eq. ^ are shown in Fig. |5fa) 
and (b), whereas the corresponding wave functions are shown in Fig. [6jc) and (d). 
Although details of the spectra differ from the simple description given in Eq. ©, the 
general structure remains similar. For the case u = and q ^ 0, which was studied in 
Ref. [8] , the spectrum consists of the discrete Landau levels within the energy window 
|e| < 2t and a quasi-continuous spectrum outside this window |e| > 2t. When a finite 
but non-resonant electric field 2u = qv = 0.2t is applied, the discrete Landau levels 
acquire dispersion inp x , as shown in Fig.^a). The quasi-continuous spectrum consists 
of a series of parallel lines corresponding to the L and R modes, well described by Eq. 
d3]). The wave functions are localized in the z direction both for the discrete Landau 
levels and the quasi-continuous spectrum, as shown in Fig. Hfc). However, when the 
resonance between the fields is achieved at u — qv — 0.2t, the spectrum changes 
drastically, as shown in Fig.[6fb). A continuous band of the R-movers corresponding to 
spectrum ((SJ) emerges, and the eigenstates become delocalized, as shown in Fig.^d). 
At the same time, the L-movers remain well-localized. 

We believe that the mechanism of resonant derealization of the wave functions 
in the crossed electric and magnetic fields, as discussed above, is responsible for the 
peaks in the differential conductance observed in the experiment. According to Eq. 
([T]), the position of the conductance peak is proportional to the magnetic field 



E = V /l = vH, 



(9) 
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where / is the effective length over which the applied voltage Vq drops. Equation 
(J9j) is consistent with the experimental Fig. [3J from which we obtain I = 1.2 nm, 
about three times longer than the interlayer spacing in graphite. This result indicates 
that decoupling of the interlayer Bernal correlation probably occurs in two or three 
intrinsic tunnel junctions. In NbSe3 mesas, it was demonstrated that the charge- 
density-wave decoupling can occur within two intrinsic tunnel junctions of the mesa 
|12) . Although our theoretical model was developed for a large number of graphene 
layers, the resonant condition (p} is expected to be valid even for a few layers. 



4 Conclusion 

The interlayer tunneling spectra dl/dV measured in graphite mesas at a strong in- 
plane magnetic field H demonstrate a peak at the voltage Vq proportional H. The 
experimental result is consistent with the presented theory, where the wave functions 
of the Dirac electrons in graphene layers delocalize in the out-of-plane direction when 
the resonant condition E = vH between the fields is achieved. 
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